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TECHNICAL NOTE 5065 


PRESENT STATUS OF INFORMATION RELATIVE TO THE PREDICTION 
OF SHOCK-INDUCED BOUNDART-LAYER SEPARATION 
By Roy H. Lange 


SUMMARY 


The present status of available information relative to the predic- 
tion of shock- induced boundary-layer separation is discussed. Experi- 
mental results showing the effects of Reynolds number and Mach number on 
the separation of both 1 ami nar and turbulent boundary layers are given 
and con5)ared with results obtained by available methods for predicting 
separation. The flow phenomena associated with separation caused by 
forward-facing steps, wedges, and Incident shock waves are discussed. 
Applications of the flat-plate data to problems of separation on spoilers, 
diffusers, and scoop Inlets are indicated for turbulent boundary layers. 


INTRODUCTION 


One of the fundamental problems that appears in the investigation of 
supersonic flow over a surface is that of the phenomena associated with 
the interaction of shock waves and boundary layers. The problem of whether 
a given shock wave will cause boundary- layer separation is one which occurs 
in all cases where a pressxrre increase is to be obtained as a result of the 
retardation of the flow. Such problems occiar, for exan^le, in the flow in 
supersonic diffusers and air Inlets and in the flow at the rear of airfoils 
and bodies. Shock- induced boundary-layer separation generally results in 
poor aerodynamic efficiency in the former case and in undesirable airfoil 
characteristics in the latter case and, therefore, this problem is of con- 
siderable practical significance. The purpose of this paper is to discuss 
the status of information, available as of July 1 , 1955 ^ relative to the 
prediction of shock- induced bovindary-layer separation. In order to study 
the fundamental featxrres of the problem, the discussion is concerned 
principally with data obtained on flat plates in two-dimensional flow. 
Exper i m en tal results showing the effects of Reynolds number and tfe,ch num- 
ber on the seperation of both 1 ami na-r and turbulent boundary layers are 
given and compared with results obtained by available methods for pre- ' 
dieting separation. The flow phenomena associated with sepaxation caused 
by f onward-facing steps, wedges, and incident shock waves are disc^issed. 
Applications of the flat-plate data to problems of separation on spoilers, 
diffusers, and scoop inlets are indicated for turbixLent boundary layers. 
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DISCUSSION OF EROEEJEM 


Erandtl has dlsciossed separation of the inconpressible boundary layer 
under the influence of a positive pressure ^gradient (refs. 1 and 2). The 
approximate methods such as those of Von Karman^ Po h l h aus en ^ and Buri for 
predicting separation were derived on the asstmoption that the boundary 
layer has time to adjust itself to a prescribed pressure distribution. 

Olhe Von Khrmidh-Pohlliausen approximation for a 1 .a m1n a r boundary layer is 


6 


dx 


^1 




-1 


( 1 ) 


and Buri’s approximation for a turbulent boundary layer is 



<ll 



( 2 ) 


boundary-layer thickness 
streamwlse pressure gradient 

free-stream dynamic pressure 
enpirical constants 

Reynolds number based on distance 6 

Experience with the xise of these approximations has shown that the occur- 
rence of separation depends chiefly upon the pressure gradient dp/dx, 
and that the turbiilent boimdary layer can withstand a much greater pres- 
sure Increase before separation than can a laminar boundary layer. 

When the influence of a shock wave on a boundary layer is considered, 
it is evident that, if the infinite free-stream pressTxre gradient which the 
shock wave represents could extend all the way to the wall, then separa- 
tion would certainly result; however, as shown by the sketch in figure 1, 
it is known that the pressure difference across the shock is spread out 
in the lower levels of the boundary layer. (See refs. 5 to J.) The work 
of Liepmann and Acker et has shown that the amount of spread of the pres- 
sure rise at the wall depends upon the state of the boundary layer, that 
is, whether the boimdary layer is laminar or turbulent (refs. 5 and 6). 
Thiis, the pressure gradient appearing at the wall boundary is fixed by the 


where 

5 

dx 


Ki,Kg 
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physical properties of the homdary layer and hy the strength of the shock 
wave. It seems logical to assume, then, that the occurrence of separation 
in this case depends principally upon the pressure rise p^ - through 

the shock wave. It was further anticipated that as the pressure rise across 
the shock was decreased there would he one shock strength below which no 
separation of the boundary layer would occur. This concept was advanced 
by Beastall and Egglnk at the British National Physical Laboratory and, 
later, a slr^jllfled dimensional analysis presented in reference 8 indi- 
cated that the critical pressure rise across the shock AP/q.^. 

just causes separation of the boundary layer shoifLd be proportional to 
the local skin-friction coefficient c^. These approximations are extended 

to the case for flat plates in terms of the Reynolds nimiber based on x. 
Thus, for a laminar boundary layer, 

^ cc cc Rg"^ cc R/^/^ (5) 


and for a turbulent boxmdary layer 

t 

^ cc Cf (^) 

It should be onphasized that the relationships given in equations (l) to 
(4 ) are only approximations . For incompressible flow more refined methods 
have been developed (refs. 9 to 12); however, the applications of these 
methods for predicting separation have met with only limited success. A 
collection of the available data for supersonic flow (ref. 8) appeared to 
bear out the predictions shown by equations (3) and (4) at the time they 
were first derived; however, since that time, more experimental data have 
come to li^t, especially for the turbulent boundary layer, which show 
that the problm must be reexamined. 


LAMINAR BOUNDARY LAYER 


The available data for shock- induced separation for ^am^■ns^r boundary 
layers on flat plates are given in figure 2, where the critical pressure 
rise across the shock divided by the free-stream dynamic pressure ^Ap 

is plotted against Reynolds number R^ on logarithmic scales. The Reynolds 


mmiber is based on the distance from the leading edge of the plate to the 
point of intersection of the shock wave and the boundary layer. The sources 
of these data are given at the top of the figure. (See refs. 3, 6, and I3. ) 
The data at Mach numbers of 1.40, 1.44, and 2.05 (refs. 3 and I 5 ) were 
obtained from tests in which shock waves of varying strength were made to 
impinge upon the boimdary layer on a flat plate. For laTnlnm- boundary 
layers which experience Interaction with incident shock waves, Liepmann 
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and Ms associates (ref. 3) liave well establislied the model for the 
resulting complex flow patterns. It has been found that the shock-wave 
patterns produce a pressure distribxrtion with the pressure rise in two 
steps, and as described in reference 8 the criticeil pressure rise is taken 
at the knee of the first step of the pressure distribution. The data at a 
Mach number of 1-95 were obtained in the separated region ahead of a 
forward-facing step. It can be seen that the available data are rather 
limited in scope and, therefore, are not conclusive; however, there are 
sane trends in the data wMch should be mentioned. For example, at free- 
stream Mach numbers of 1.93 and 2.O5 the Reynolds number effect on 

the critical pressure coefficient appears to follow the inverse square 
root of the Reynolds number as denoted by the dashed lines in figure 2. 
Except for the data at Mach numbers of 1.40 and 1.44, the critical pres- 
sure coefficient also decreases with increasing Mach number. These trends 
of Reynolds number emd Mach number agree with the predictions of equa- 
tion (3); however, the magMtude of the Mach number effect shewn, espe- 
cially between Mach numbers of 1 and 2, is much greater than that which 
would be predicted by reference 8. Recent data obtained at the Ames 
Laboratory in the separated region ahead of a forward-facing step show 
an increase in with increase in Reynolds number wMch is 

opposite to the trend obtained at a Mach number of 1.93* The forward- 
facing-step data shoim in figure 2 should be ignored, therefore, until 
more systematic data are available. Stewart son (ref. l4) has made a 
detailed analysis of the interaction process wMch leads to the inference 
that the dimensionless pressure rise required to produce laminar separa- 
tion would be proportional to -2/5. 

Also shown in figure 2 is a curve 

wMch traces the criterion of separation advanced by Pabst (ref. I3) in 
a recent Argentine paper; however, tMs criterion cannot account for the 
Mach number effect and does not correlate with any of the experimental 
data shown. 


TURBULENT BOUNDARY LAYER 

O 


Investigations of shock — ^boundary-layer Interaction for the tiarbulent 
boundary layer have shown that a given shock wave may or may not separate 
the boundary layer. Data are now available from a number of sources in 
wMch turbulent-bo\mdary-layer sepairation has been investigated by three 
methods; (l) the forward-faclng-step technique, (2) the wedge techMque, 
and (3) the incident-shock techMque. 

In order to remove all doubt as to whether the turbulent boundary 
layer has been separated, several investigators have forced separation 
by means of a forward-facing step moxnited on a flat plate (see refs. 8, 
16, and 17) . "lypical data for tMs type of configuration are given in 
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figure 5 'wMcli shows the pressure distribution along the surface and 
(to the same scale) a sketch of the flow field in the interaction region 
as determined from shadowgraphs. These data were obtained in a blowdown 
jet of the Langley gas dynamics laboratory at a Mach nuniber of 3* 05* The 
model and supersonic nozzle were the same as those described in detail in 
reference 8 except that the model was equipped with static-pressure ori- 
fices. The flow diagram at the top of the figure shows that a wedge- 
shaped separation region is fonaed ahead of the step and is bounded on 
its upstream edge by the shock wave. The direction of the circulatory 
flow within the separated region is shown by the arrows. The pressure 
coefficients on the plate first reach a maxltnum value, noted herein as 
the first peak, at a point about halfway between the location of the 
shock wave and the location of the step. This distance is roughly the 
equivalent of 8 boundary-layer thicknesses or 133 momentum thicknesses, 
on the assumption of a 1/7-power velocity distribution In the boundary 
layer just ahead of the shock. The pressures then dip sli^tly behind 
the first peak and subsequently rise shaiply; this result shows the large 
influence of the circulatory flow. 

Also pertinent to the discussion of the flow in the separated region 
are the pressure coefficients measured along the front vertical face of 
the step given in figure 4. The data for Reynplds numbers ranging from 

12 X 10^ to 32 X 10^ were obtained at = 3*05 at the Langley gas 

dynamics laboratory. The pressure orifices were located at the base of 
the step and at two other vertical locations above the surface of the 
plate as denoted by z/h. The data at Mi = 3*03 show no significant 
Reynolds nimiber effect on the pressure coefficients. The results show 
that there is one stagnation point at the foot of the step and one near 
the top of the step, and calculations based on the data at Mi = 3*03 

and utilizing the inconpressible Bemoul li equation show that the veloc- 
ity downward along the vertical face is about one-fourth the free-stream 
velocity; whereas the velocity along the plate in a direction opposite 
to the main flow is about one-third the free-stream velocity. Thus, the 
separated region cannot be treated as a dead-air space as is commonly 
assumed. The results also show that a considerable error would result 
if the pressiires on the front face of the step were assumed to be the 
same as those obtained on the plate surface ahead of the step in the 
separated region. The first peak pressure coefficient obtained ahead 
of the step is shown by the dashed line at a Mach nimiber of 3 •05* lb 
is clear then, from the res-ults given in figures 3 an^ that the first 
peak pressure coefficient is obtained as a result of the mutual effects 
of the shock on the boundary layer and of the circulatory flow in the 
separated region and should not be interpreted as the value of the pres- 
sure rise across the Tnjnimum strength of shock wave which just causes 
separation of the boundary layer. These results have been obtained for 
cases where the step height is about 3 times the local boundary-layer 
thickness and may be changed somewhat for cases where the step hei^t 
is very large compared with 5. 
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A summary of tlie available data obtained from tlie Tise of the step 
technique for forcing boundary-layer separation is given in figure 5 
■which shovs taken at the first peak: plotted aga±ost Reynolds 

number on a logarithmic scale. The Reynolds number is based on the dis- 
tance from the leading edge of the plate to the point of intersection of 
the shock wave -with the boundary layer. All the data -were obtained from 
pressure distributions (see ref. 1T)> and the sources are given at the 
top of the figure. The Mach number range of the data is from 1.55 shown 
by the long string of points at the top of the data to 5*65 shown by the 
lowest data points. The pres sure -distribution data at M]_ = 5*05 given 
by the circles are new da"ta "which have not been published. The data 
given in reference 8 (MCA TR 2770) for = 3*03 represented by the 


dashed line which varies as were obtained by measuring shock 

angles close to the point of intersection of the shock "wave and the bound- 
ary layer, 'where, as shown previously, the pressures on the plate are 
changing rapidlyj therefore this method for obtaining pressxu:e coeffi- 
cients is too crude and the data should be ignored. It is apparent from 
the data that the Reynolds number effect on "the value of ('^/‘3l)ls-t peak 

is very si 1 gbt . On "the basis that there is no Reynolds number effect, 
figure 6 has been prepared to show the decrease in /q 2 j ^ 

' ' • XSu p0£LiC 

with increase in free-stream Mach nunber for Mach numbers between 1.55 
and 3«65. All the data from the previous figure have been inclioded in 
this plot, and the -vertical line connecting two of "the synbols shows the 
extent of the Reynolds number effect obtained. 


The second technique for producing "turbulent boundary-layer separa- 
tion is the use of wedges of different angles mounted on flat plates, 
and a limited amount of da"ta is available. (See refs. I 7 and I 8 . ) This 
configuration is analogous "1x3 the deflection of a flap or a control sur- 
face. [IJypical data ob"tained at a Mach nunber of 3«03 are gi"ven in fig- 
ure 7 "Which shows the pressuire distribution along the plate and on the 
wedge and abo"ve it a sketch of the flow phenomena as determined by shadow- 
graphs. A double scale is given e0.ong the abscissa of the pressure dis- 
tribution - one which gi"ves x in inches measured from the leading edge 
of the "wedge and one which gives a measure of the boundary-layer thick- 
ness x/fe. As shown in "the flow picture, the separation in the comer 
produced by this particuilar wedge angle results in a weak shock wave , 
which projects ahead of the main shock, and an inflection point is 
obtained in the pressure distribution on the surface. Downstream of 
this point the press"ure coefficient continues to rise and levels off at 
a value somewhat less than "that calculated from oblique-shock "theory for 
"this "wedge angle in "the absence of a boundary layer. In general, the 
limi"ted available data at a given Mach number show "that, for wedge angles 
greater than a certain value, the pressure distribution has an inflec- 
tion point similar "to that sho"wn in figure 7; moreover, "the value of /^p/q^ 

measured at the inflection point remains almost constant wi"th fur"ther 
increases in wedge angle. The data at M]_ = 3»03 also show that the 
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value of ^ obtained at the inflection point is essentially constant 
for Reynolds numbers ranging from 12 x 10^ to 32 X 10^. 

Results are available from tests utilizing the thlord technique 
in which shock waves of varying strength are made to in^jlnge iq)on the 
boundary layer on a flat plate. (See refs. 19 and 20.) A sketch of 
the resulting flow phenomena when the incident shock-wave strength is 
sufficient to cause local separation of the boundary layer is constructed 
from the work of Bardsley and Mair (ref. 2l) : 



As shown in the sketch, the pressure rise across the incident shock wave 
is propagated i: 5 )stream in the boundary layer and results in a rapid 
thickening of the boundary layer and in liie formation of a fan of compres- 
sion waves just ahead of the point of intersection of the Incident shock 
wave and the boundary layer. The incident shock wave is reflected from 
the boundary layer as an expansion (as at a constant-pressure sttrface). 
This shock-expansion process turns the flow toward the plate and the rate 
of thickening of the boundary layer is decreased. Downstream of this 
point a second set of con^iression waves is formed which ultimately returns 
the flow parallel to the plate. The pressure alorg the plate increases 
rapidly tip to the point where the flow separates from the plate, after 
which the pressure gradient decreases until the second compression fan 
Increases the pressirre to the tiltimate value. Thus the pressure distri- 
bution has an inflection point similar to that obtained in the wedge 
tests, and the pressure rise for separation is taken at the inflection 
point. 
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It Is interesting to note that the shock-wave — honndary-layer Inter- 
action process just described for the turbulent boundary layer is simi- 
lar to that given by Liepmann for separation of the laminar boundary 
layer; however, the magnitudes of the upstream influence of the pressvire 
rise across the incident shock waves and of the press-ure coefficients 
for separation are considerably different for the laminar and turbulent 
boundary layers. 

In figure 8 the values of for separation of the turbulent 

boundary layer are plotted against Mach number, vhere the inflection- 
point pressure coefficients obtained in the wedge tests are given by 
the open symbols and the inflection-point pressure coefficients obtained 
by the incident-shock technique are given by the solid symbols. The 
experimental data point shovn at M = 2 for the incident-shock tech- 
nique was obtained at the British National Physical Laboratory by Gadd 
and Holder as imported by Li^thill in reference 22. Also shown in this 
figure is the curve representing the data obtained by the forward-faclng- 
step technique. The data given in this figure, therefore, constitute 
all Information available at present on turbulent-boundary-layer separa- 
tion. The spread in Z?)/q]_ obtained at = I .80 in the wedge tests 

represents a Reynolds number effect, althou^, as mentioned previously, 
no such Reynolds number effect was obtained at = 3*03» Althou^ the 
available data are rather limited in scope, the results show that the 
inflection-point pressure coefficients obtained from both techniques gen- 
erally have the same range of values with increasing Mach number and that 
on the average these values are about 20 percent lower than those obtained 
by using the step technique. The application of these data for predicting 
separation should, therefore, be limited to these particxiLar configurations, 
at least for the present. For example, the data from the incident-shock 
technique represent conditions of local separation of the flow and, because 
the experiments are performed on flat plates, the flow reattaches down- 
stream of the separation point. This reattachment may be changed some- 
what for conditions where a back pressure exists - for example, for con- 
ditions near the trailing edge of an airfoil. Also, fli^t data for a 
wing in transonic flow Indicate that the value of isg for separation 

is predicted more accurately by the step data if extrapolated to the 
lower st^rsonlc Mach numbers obtained in the fli^t tests (ref. 23)* 

These data are useful, then, in providing a first approximation to the 
pressure coefficient for which separation is likely to be encountered. 


CONCLUSIONS 


In conclusion, the present status of information relative to the 
prediction of shock-induced boundary-layer separation indicates that, 
althou^ no -universal value of pressure-rise coefficient which causes 
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Incipient separation of the houndary layer has heen found, there is a 
feiirly narrow hand of pressure coefficients from ■which predictions of 
turbulent separation can he made with an accuracy prohahly sufficient 
for engineering purposes. On the basis of these results the following 
tentative conclusions are given: 

1. Ihe da'ta obtained 'with fonfard-facing steps, -wedges, and inci- 
dent shock waves indicate -that there is a dependency of -the pressure 
coefficient for separation on Reynolds number for the laminar boundary 
layer but little, if any, dependency on Reynolds number for the turbu- 
lent boundary layer. Qiiere is a dependency of this pressure coefficient 
on Mach nuniber for both laminar and 'turbulent boundary layers. 

2. For application to. a spoiler, the available da'ta obtained by the 
forward-facing-step technique permit calculations of the loading on the 
surface ahead of the spoiler, the pressure on the front face of the 
spoiler, and the separation point ahead of the spoiler for a Mach num- 
ber range from 1-55 "to 3*65 for the turbulent boundary layer. 

3 . For application to supersonic diffusers or scoop inlets, the 
available data from incident-shock-wave tests pro'vlde a first approxi- 
mation to the mi-n-tTmim strength of shock -which will separa-be the turbu- 
lent boundary layer for Mach numbers be-bween 2 and 3* 

If. From the da-ba available from the wedge -bests, a first approxi- 
mation to -bhe pressiire coefficient for -which separation becomes apprecia- 
ble as a result of flap deflection can be made for a surface -with a -turbu- 
lent boimdary layer for Mach numbers between 1.75 3* 03* 

5 . Caution sho-uld be exercised in attempting to predict the separa- 
tion or loading on configurations which differ considerably from -those 
for -which experimental data are available. 


Langley Aeronautical Labora-tory, 

National Advisory Commi-ttee for Aeronautics, 
Langley Field, Va. , Septeiriber 1, 1955* 
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Figure 1 .- Shock wave entering a boundary layer. 
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Figure 2 .- Critical pressure coefficients for laminar boundary layers. 
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Figure 5 
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Figure 7 *- Separation of turbulent boundary layer by wedge. Data 
obtained in L an gley gas dynamics laboratory. 
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Figure 8 .- Separation pressure coefficients for turbulent boundary layers 
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